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Abstract 

This paper addresses the stability and queueing delay of Space Division Multiple Access (SDMA) 
systems with bursty traffic, where zero-forcing beamforming enables simultaneous transmission to 
multiple mobiles. Computing beamforming vectors relies on quantized channel state information (CSI) 
feedback (limited feedback) from mobiles. Define the stability region for SDMA as the set of multiuser 
packet-arrival rates for which the steady-state queue lengths are finite. Given perfect CSI feedback and 
equal power allocation over scheduled queues, the stability region is proved to be a convex polytope 
having the derived vertices. For any set of arrival rates in the stability region, multiuser queues are 
shown to be stabilized by a joint queue-and-beamforming control policy that maximizes the departure- 
rate-weighted sum of queue lengths. The stability region for limited feedback is found to be the perfect- 
CSI region multiplied by one minus a small factor. The required number of feedback bits per mobile is 
proved to scale logarithmically with the inverse of the above factor as well as linearly with the number 
of transmit antennas minus one. The effects of limited feedback on queueing delay are also quantified. 
For Poisson arrival processes, CSI quantization errors are shown to multiply average queueing delay 
by a factor M > 1. For given AI — > 1, the number of feedback bits per mobile B is proved to be 
0{— log2(l — )). For general arrival processes, CSI errors are found to increase Kingman's bound 
on the tail probability of the instantaneous delay by one plus a small factor rj. For given 77 0, B is 
proved to be 0( — logj ij). 



I. Introduction 

In this paper, we consider a Space Division Multiple Access (SDMA) system where a multi- 
antenna base station transmits simultaneously to multiple single- antenna mobiles. Given feedback 
of channel state information (CSI) from mobiles, data packets transmitted using SDMA are 
decoupled by zero-forcing beamforming at the base station. To design such a system, it is 
important to understand the maximum packet arrival rates the system can support without 
becoming unstable, and the effects of the practical finite-rate CSI feedback (limited feedback) on 
system stability and queueing delay. We address these issues by deriving the stability region for 
zero-forcing SDMA, defined as the set of packet arrival rates for which finite steady-state queue 
lengths are feasible. We also characterize the required amount of limited feedback overhead for 
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constraining the degradation of system stability and delay performance due to limited feedback. 
These results provide insight into designing admission and queue control as well as limited 
feedback in multi-antenna SDMA systems. 

A. Prior Work and Motivation 

SDMA for multi-antenna downlink systems has emerged as a key technology for enabling 
high-rate wireless access [1]. The key feature of SDMA is the support of multiuser data streams 
by exploiting the spatial degrees of freedom. This feature can be realized by pre-cancellation of 
multiuser interference using the optimal dirty paper coding (DPC) [2], [3]. This technique not 
only has high complexity but also requires non-causal and perfect CSI at the transmitter (CSIT). 
For these reasons, multiuser beamforming such as zero-forcing beamforming has become a 
popular alternative solution by providing low complexity and close-to-optimal performance [4]- 
[7]. Furthermore, multiuser beamforming admits efficient scheduling [8]-[10]. 

In practice, feedback CSI for multiuser beamforming must be quantized given finite-rate 
feedback channels [11]. Limited feedback introduces interference between multiuser data streams 
even if zero-forcing beamforming is applied [5], [9]. The resultant throughput loss can be 
controlled by adjusting the number of CSI feedback bits. In particular, for zero-forcing beam- 
forming, this number must scale with the number of transmit antennas and the signal-to-noise 
ratio (SNR) to constrain the throughput loss [5]. If the user pool is large, the effects of CSIT 
quantization errors can be alleviated by exploiting multiuser diversity [6], [9], [10], [12]. The 
above works on SDMA and limited feedback assume mobiles with infinite back-logged data. 
However, it is important to address the issue of data burstiness for SDMA supporting delay- 
sensitive applications. 

Active research has been conducted on multiuser systems with bursty traffic by integrating 
information and queueing theory. In [13], the delay and throughput of multiple-access channels 
with Poisson packet arrivals are studied as a processor- sharing problem, where data packets 
are decoded by treating multiuser interference as noise. In [14], [15], the more sophisticated 
successive decoding [16] is applied in multiple access channels with bursty data. These studies 
show the delay and stability optimality of the queue control policy that allocates the highest data 
rates to mobiles with the longest queues, aligned with the well-known principle of max weight 
match [17]. Given successive decoding, the stability regions for both the multiple-access and the 
broadcast channels are found to be identical to their corresponding information capacity regions 
[18]. A general broadcast system of multiple queues and multiple servers is studied in [19], 
where the max-weight-match policy for resource and server allocation is shown to be stability 
optimal. These works assume perfect CSIT available at the base station. The issues of stability 
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and queueing delay for SDMA with limited feedback are challenging and still unexplored. 

Existing work on multi-antenna SDMA with bursty data focuses on enlarging the stability 
region by exploiting multiuser diversity with perfect CSIT [20], [21]. Based on the max- weight- 
match principle, these scheduling algorithms select mobiles with long queue lengths and whose 
channel vectors are nearly-orthogonal with each other. The orthogonality criterion reduces 
multiuser interference [21] or signal power loss due to interference avoidance [20]. However, 
the existing algorithms are ineffective for systems with small numbers of mobiles, or those 
supporting delay-sensitive applications where scheduling based on mobiles' CSI is infeasible 
[5]. For these systems, we propose the approach of integrating queue control and adaptive 
beamforming. Thereby the downlink stability region is enlarged by dynamic control of the 
number of active queues and spatial streams. 

B. Contributions and Organization 

We consider a narrow-band SDMA system where a multi-antenna base station serves multiple 
single-antenna mobiles using zero-forcing beamforming with limited feedback. The number 
of mobiles is equal to that of transmit antennas. The channel coefficients are i.i.d. CJ\f{0, 1) 
random variables. The data packets for different mobiles are stored in separate buffers of infinite 
capacity at the base station, forming multiple queues. In each time slot, a subset of queues are 
scheduled and corresponding data packets are transmitted using zero-forcing beamforming and 
automatic repeat request (ARQ) [22]. Packet transmission is successful if the received signal- 
to-interference-and-noise (SINR) exceeds a given threshold. 

We investigate the stability and delay performance of the above system. Our main contributions 
are summarized as follows. Given perfect CSIT and equal power allocation over scheduled 
queues, the SDMA stability region is proved to be a convex polytope having the derived vertices. 
Each vertex corresponds to scheduling a unique subset of queues for SDMA. Queues with arrival 
rates in the stability region are shown to be stabilized by the queue-and-beamforming control 
policy that maximizes the departure-rate- weighted sum of queue lengths. In addition, we obtain 
the stability region for flexible power control. Next, we analyze the feedback requirements 
for zero-forcing forming with limited feedback. The required number of feedback bits per 
user B is derived under the constraint that the stability region is no smaller than the perfect- 
CSI counterpart scaled by one minus a small factor. Specifically, B is shown to increase 
logarithmically with this factor as well as linearly with the number of transmit antennas minus 
one. The feedback requirements are also analyzed using the criterion of queueing delay. For 
Poisson arrival processes, B is obtained for bounding the average-delay ratio between limited 
and perfect CSI feedback. In particular, B is shown to have the same order as the logarithm of 
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Fig. 1. Multi-antenna SDMA system 

the above ratio if it is close to one. For general arrival processes, the complementary cumulative 
distribution function (CCDF) of instantaneous delay is upper bounded by Kingman's bound [23]. 
Using perturbation theory [24], B is derived given a constraint on the multiplier of the CCDF 
bound caused by CSIT inaccuracy. 

The remainder of this paper is organized as follows. The model of the SDMA system is 
described in Section [III In Section [nil assuming perfect CSI feedback, we discuss the stability 
region for SDMA and the joint queue-and-beamforming control policy for stabilizing the SDMA 
system. The CSI feedback requirements for constraining the stability and delay performance loss 
due to imperfect CSI are analyzed in Section |IVl Numerical results are presented in Section |Vl 
followed by concluding remarks in Section |VIl 

II. System Description 

We consider a narrow-band SDMA system where a base station with L antennas serves L 
single- antenna mobiles as illustrated Fig. [T] Following [5], [25], [26], channel coefficients are 
modeled as i.i.d. CA/'(0, 1) random variables, corresponding to a rich scattering environment. 
Channel variation over time is modeled as block fading, namely that each realization of channel 
coefficients is fixed within one time slot and different realizations are independent [27]. Without 
loss of generality, each time slot is assumed to span one time unit. Moreover, the fading process 
is assumed stationary and ergodic. The vector channel for the £th user is represented by a L x 1 
random vector h.(^{t) where t is the slot index. 

A. Data Source and Transmission Protocol 

The data for each user is generated in packet. For simplicity, the packet length is assumed 
fixed and equal to one time unit. Packets for different mobiles arrive asynchronously at separate 
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buffers with infinite capacity, fomiing L queues at the base station. Multiuser packet arrival rates 
may be unequal, modeling multimedia applications. At the beginning of each slot, packets from 
scheduled queues are transmitted using zero-forcing beamforming (see Section III-BI) . Packet 
transmission uses the following ARQ protocol [22]. A packet is received successfully if the 
SINR at the intended mobile exceeds a given threshold 6. Otherwise the packet is retransmitted 
in the next scheduled time slot. Each mobile reports to the base station the status of each received 
packet via perfect ACK/NACK feedback. 

B. Zero-Forcing Beamforming 

Consider a set of K nonempty queues scheduled for transmission in an arbitrary time slot, 
which are specified by the index set A. For simplicity, we assume equal allocation of the total 
transmission power P over scheduled queues, and the effect of relaxing this assumption on 
system stability is discussed in Section IIII-CI The K transmitted packets are separated using 
beamforming vectors represented by {ff : £ G A]. Let Xi := CJ\f{0, 1) and Zc represent the 
transmitted and received data symbols for the ith user, respectively, where the symbol index is 
omitted for brevity. We can write Zi as 

where Ni := CJ\f{0, 1) denotes a sample of the additive white Gaussian noise process, and 
7 := ^ is the signal-to-noise ratio (SNR). Note that the summation term in ([T]) represents multi- 
user interference. Let he denote the feedback CSI generated from h^. The base station computes 
beamforming vectors using {h^ : C, G A} and the zero-forcing method [5]. To be precise, is 
chosen to be orthogonal to the channel subspace He spanned by the vectors {hm \ m i,m E A} 
[5], [9]; under this constraint, maximizes the received SNR, namely fe = argmaxfgT^c |f/h£p 
where He denotes the null space of H [28]. For perfect CSIT, these zero-forcing constraint nulls 
multi-user interference (the summation term in ([1]) is equal to zero) for the case of perfect CSIT 
{he = he) or otherwise residual interference exists. Thus the SNR/SINR can be written from ([T]) 
as 

' S N = 7 1 f j 1 2 perfect CSI feedback 

< 7|fJh.P (2) 
SINR^ = n e i\ limited feedback. 

Zero-forcing beamforming with perfect CSIT and limited feedback are considered in Section UlI] 
and Section |lVl respectively. 
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C. Limited Feedback 

In practice, zero-forcing beamforming relies on limited feedback from mobiles. Specifically, 
each mobile estimates the corresponding vector channel using pilot signals broadcast by the base 
station. Then all mobiles send back CSI via finite-rate feedback channels to the base station for 
enabling multiuser beamforming. We assume perfect CSI estimation at mobiles and zero delay 
and error for the feedback channel following the literature (see e.g. [5], [29]). Thereby the 
CSIT inaccuracy is contributed entirely by CSI quantization. To simplify analysis, we adopt the 
following quantized CSI model from [9], [30]. Let a unitary vector represent the output of 
quantizing h^. Define the quantization error as = 1 — |hjh^p/||h^p. The model in [9], [30] 
approximates a Voronoi cell [31] as a sphere cap, or more precisely is isotropic in the set 
(Voronoi cell) V = |s G : ||s|| = 1, < 1 - 2"A|. As a resuk, the CDF of is [9] 

Pr(e^ < a) = 2^a^"\ < a < 2"^. (3) 

The above model is observed in [9], [30] to be accurate for the practical codebook-based CSI 
quantization for beamforming [32]. 

D. Joint Queue-and-Beamforming Control 

In each slot, the controller (see Fig. [T]) schedules a subset of nonempty queues for SDMA 
transmission. Furthermore, the controller informs the beamformer to create the matching spatial 
streams for transmitting packets from scheduled queues. For simplicity, we consider equal 
transmission power for scheduled queues. This assumption is relaxed in Section IIII-CI Under 
this assumption, the controller's decision in the tth slot can be represented by an L x 1 indicator 
vector m(t) G V := {0, 1}^, where the £th component of m(t) is equal to one if the £th queue 
is scheduled or otherwise equal to 0. Let Qi{t) G IN denote the £th queue length and define the 
vector q(t) := [Qi{t),Q2{t), ■ ■ ■ .Qiit)]. The policy for joint queue-and-beamforming control, 
referred to hereafter simply as the control policy, is defined as a function tt : IN^ — > V, and 
thus m(t) = 7r(q(t)). A stationary policy refers to one for which Pr(m = v) for any V is time 
invariant. The control policy is further discussed in Section IIII-B[ 

III. Stability for Perfect CSI Feedback 

In this section, assuming perfect CSI feedback, we analyze the stability region for zero-forcing 
SDMA and discuss joint queue-and-beamforming control for stabilizing queues. Consider equal 
power allocation over scheduled queues. The stability region is shown to be a convex poly tope 
with derived vertices. For any arrival-rate vector in this region, queue stability is shown to be 
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achieved by the control policy that maximizes the departure-rate-weighted sum of queue lengths. 
Finally, we discuss the stability region for unequal power allocation. 

A. Stability Region 

The departure rates for the SDMA system are derived. Then its stability region is defined and 
shown to be a convex poly tope. 

Several useful definitions are provided. Let /x(m) denote the L x 1 departure-rate vector 
conditioned on the control decision m (see Section III-DI) . The ith component of /x(m) gives 
the departure rate of the £th queue and is defined as [/x(m)]^ := Pr(SNR^ > 9 \ m) where SNR^ 
is in (l2l). Let denote the ith of the 2^ element in the control decision space V. A stationary 
policy TT can be specified by a probability vector = [pi,P2, ■ ■ ■ ,^2^] where = Pr(m = v^) 
and IIPttIIi = 1 (see Section HLDI) . § Let /i(p7r) denote the departure-rate vector conditioned on 
TT. It follows that /i(p^) = Y.'eli m(v<?)P£. 

The conditional departure-rate vector //(Ptt) is derived as follows. Let K := ||m||i represent 
the number of scheduled queues. Zero-forcing beamforming for fixed K < L is considered 
in [28] and called partial zero-forcing beamforming. From [28, Lemma 2], given K = k, the 
effective channel power for a scheduled queue, namely |f/h^P in follows the chi-squared 
distribution with {L — k + 1) complex degrees of freedom, denoted as — A; + 1). Therefore 
the departure-rates for all scheduled queues are identical and given by d{k) in (|5]). Note that 
d{k) decreases with increasing k and vice versa, where k also specifies the spatial multiplexing 
gain. It follows that /x(m) can be written as /x(m) = (i(||m||i)m. Using this result, /i(p7r) is 
obtained as shown in the following lemma. 

Lemma 1: Given perfect CSIT and the stationary control policy vr, the conditional departure- 
rate vector /i(p,r) is given as 



KPtt) = 2^ Pidi\\ve\\i)vi (4) 



where d{k) given below is the departure rate for each of total k scheduled queues 

, , r(L-k + i,!^) 

= r(L-.H-i) ■ <5> 

Using these results, we now define the stability region for zero-forcing SDMA. Adopting the 
definition in [33], the ith queue at the base station is stable if the queue length Qe{t) (in packet) 

'll^lli ~ denotes the Li norm of the vector a. 

^r(a) ~ JJ" t"^^e^'dt and r{a,b) := JJ^ f^^^e^^dt denote the gamma and the upper incomplete gamma functions with 
a complex degrees of freedom, respectively. 
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satisfies 




lim Pr(Q^(t) < q) = F{q) and lim F{q) = 1 



(6) 



wliere F(-) denotes a CDF. Let and fig denote the arrival and departure rates for the ith queue, 
respectively. It follows from Loynes' theorem [34] that the £th queue is stable if < ftg, but the 
stability condition for the boundary point A^ = fig is uncertain. To simplify our discussion, we 
assume stability at A^ = jl£ as in [35]. Thereby we can define the stability region for the policy 
TT as the closure of arrival-rate vectors for which all queues are stabilized by tt. To be precise, 
given 71 and the corresponding probability vector p^^, the stability region ^(p^) is defined as 



where ^ represents component-wise inequality and fJ-iPn) is in Lemma [T] The stability region C 
for SDMA is readily defined below as the union of the stability regions for all feasible policies 



where ^(p^) is given in ([7]). 

Define the set TZ := {/x(v) : v G V} = {rf(||v||i)v : v G V} that groups the departure-rate 
vectors for all control decisions in the space V. The following lemma states that the stability 
region C is a convex polytope whose vertices belong to TZ. 

Lemma 2: Given perfect CSIT, the stability region for SDMA is C = c6 (TZ). § 
Proof: See Appendix [2l □ 

Lemma [2] shows that C is a convex polytope with a subset of points in TZ being the vertices. 
Other points are not extreme points of C. They lie either on the surface of the polytope C or 
inside it. 

We are now ready to specify the vertices of the stability region C. To this end, some useful 
notation is introduced. Conditioned on K = k, the control decision m belongs to the set 
Vfc := {v G V : llvp = k}. The corresponding set of departure-rate vectors is defined as TZk : = 
{d{k)v : V G V^} where d{k) is given in ©. Note that Uo<fc<L^fc = ^' Uo<A:<L^fe = 
Vo = 7^0 = {0}. Define the index set 



AiPn) := {A G I A ^ /x(P-)} 



(7) 




(8) 



Pir||l=l 




(9) 



The main result of this section is given in the following theorem. 



CO denotes the closed convex-hull operation [36]. 
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Fig. 2. Stability region for L = 3 and (a) low SNR: P/O = 0.5 or (b) high SNR: P/O = 10 

Theorem 1: For perfect CSIT, the stability region C is a convex polytope in whose vertices 
(extreme points) are given in the following set 



Proof: Essentially, we prove this theorem by showing that ext C contains all the exposed points 



Note that ext C contains the origin 0. Theorem \T\ suggests the following simple procedure for 
constructing C, where -Rmax denotes the maximum downlink packet rate. 

1) Initialization: k = 0, ext C = {0}, and -Rmax = 0; 

2) Let k = k + 1. If kd{k) > /?max, ^max = kd{k) and ext C = ext C U {7^fc}; 

3) If A; < L, repeat 2); otherwise go to 4) 

4) Construct a polytope C using the points in ext C as vertices. 

Step 2) implies that scheduling k queues enlarges the stability region only if it increases the sum- 
packet rate with respect to scheduling fewer queues. Note that the number of scheduled queues 
achieving -Rmax varies with the SNR and need not be equal to the maximum L. In other words, 
the shape of the stability region C varies with SNR as illustrated by the following example. 

Example 1: The stability region for -L = 3 is plotted in Fig. |2a) for P/6 = 0.5 and Fig. |2tb) 
for P/6 = 10, corresponding to a low and a high SNRs, respectively. As observed from Fig. [2l 
ext C reduces to ext C = T^i U {0} for the low SNR and has the maximum size, namely 
ext C = TZ, at the high SNR. This observation is formally stated in the following proposition. 




(10) 



of C. □ Details are presented in Appendix |Bl 



□ 



''a point g G C is exposed if there exists a hyper-plane V supporting C such that V HC = g [36]. 
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Proposition 1: There exist r > such that ext C = TZiU {0} \/ P < t, and Pq > such that 
ext C = 7^ V P > Po- 

Proof: See Appendix O □ 

The above proposition suggests that at low SNRs the control policy based on time-division- 
multiple-access (TDMA) is stability optimal, namely stabling queues for any arrival-rate vector 
in the stability region; at high SNRs, the stability-optimal policy should schedule multiple queues 
for SDMA. Discussion on the control policy is presented in the next section. 

B. Joint Queue and Beamforming Control 

In this section, we discuss stability optimal policies for joint queue-and-beamforming control 
(see Section ini). 

If the arrival rates of the queues are known, the queue stability can be achieved by pre-defined 
time-sharing between scheduling different subsets of queues. An arrival rate vector A G C can be 
written as a convex combination of the vertices of C since it is a convex polytope. Specifically, 
A = Yl^=i'^n^n, where u„ is the nth element of ext C in (fTOl) . T := |ext C\u and {a„} are 
positive and satisfy the constraint Yln=o = 1- We can find at least one point A' on the boundary 
of C such that A' y A. Write A' = Yln=i ^n'^n with 6„ > V n and Yln^n = 1- Recall from 
earlier discussion that we can write u„ = (i(||v„||i)v„ where v„ G {0, 1}^. Therefore, given A, 
queues can be stabilized by any time-sharing policy where the fraction of time when m(t) = v„ 
is bn- Alternatively, a randomized policy can be applied such that Pr(m = v„) = 6„ [17]. 

In practice, the packet arrival rates are usually unknown. For this case, the queues can be 
stabilized using a policy that controls queue and beamforming based on queue lengths rather 
than packet arrival rates. Using the principle of maximum weight matching [17], the following 
proposed control policy maximizes the departure-rate-weighted sum of queue lengths 

TT : m.(t) = argmax(i(||v||i)v"^q()f:). (11) 

vGV 

The following proposition states the optimality of the above policy. 

Proposition 2: The control policy in (fTTI) stabilizes all L queues if the arrival rate vector lies 
in the stability region C in Theorem [T] 

Proving Proposition |2] uses the technique of Lyapunov drift [17] and follows a similar procedure 
as that for [19, Theorem 3]. The proof is omitted for brevity. 



^\X\ represents the cardinality of the set X. 
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Fig. 3. Effects of power control and limited feedback on the stability region for L — 2 

C. Stability Region with Power Control 

Equal power allocation over scheduled queues is assumed in the preceding sections. In this 
section, we relax this assumption and derive the stability region with power control. Compared 
with equal power allocation, power control provides an additional degree of freedom for enlarging 
the stability region as illustrated in Fig. [3l Some notation is introduced as follows, where the slot 
index t is omitted for brevity. Denote the power vector as p = [Pi, P2, - ■ ■ , -Pl]^ where J2e Pi = 
P and Pi gives the transmission power allocated to the £th queue. By extending Lemma [B the 
departure rate for the scheduled queue with the index I is (i(P^, m) := F — ||m||i + 1, -^j x 
[m]^. We define the departure rate vector /x(p, v) with the i component being (i(P£,m). The 
stability region for SDMA with power control is shown in the following proposition. The proof 
follows a similar procedure as that for Lemma [21 and is thus omitted. 

Proposition 3: The stability region for a zero-forcing SDMA system with power control is @ 



Unlike C, Cp is not a convex polytope. Moreover, C G Cp. The stability-optimal control policy 
in (fTTI) can be easily extended to include power control. 



The stability analysis in the preceding section assumes perfect CSIT. In this section, we 
consider limited feedback and analyze the feedback overhead required for constraining the 
performance degradation due to imperfect CSIT. The required number of feedback bits per 




(12) 



IV. CSI Feedback Requirements 



*For the vector v, diag(v) denotes a matrix with diagonal elements taken from v and other elements being zeros. 



Huang and Lau: Stability and Delay of Zero-Forcing SDMA with Limited Feedback 



12 



mobile B is obtained for achieving the stability region close to the perfect-CSI counterpart. We 
also derive B required for approaching optimal delay performance. 

A. Stability Region for SDMA with Limited Feedback 

CSI quantization errors not only reduce received signal power but also cause mutual interfer- 
ence between packets transmitted using SDMA. Consequently, the departure rate of each queue 
decreases, shrinking the stability region as illustrated in Fig. [3l In this section, we derive B 
for bounding the difference in stability region between the cases of perfect CSIT and limited 
feedback. 

The power distributions of the received signals and multiuser interference are given in the 
following two lemmas. Lemma [3] provides the distribution of the signal term If/h^p of the 
SINK in ©. 

Lemma 3: The random variable |f|h^P follows the exponential distribution with unit mean 
for the number of scheduled queues K = L. For K < L, |f|h£p := [eiZ — peQ), where is the 
CSI error, ee < pe < e^, and the random variables Z and Q follow independent x^{L — K + 1) 
and x^i^ ~ 1) distributions, respectively. 

Proof: See Appendix |D] □ 

This lemma generalizes [28, Lemma 2] for perfect CSIT to include the case of limited feedback. 
In other words, substituting = into Lemma [3] gives [28, Lemma 2]. 

The results in Lemma [3] are explained as follows. For full spatial multiplexing K = L, all 
spatial degrees of freedom contribute to suppressing interference between queues. Consequently, 
the beamforming vector for each queue is chosen independent of the corresponding vector 
channel [5]. That is to say, CSI is not used for strengthening received signal power and thus it 
is independent of the CSI quantization error as observed from Lemma [3l For K < L, a fraction 
of the spatial degrees of freedom (L — K + 1) provide diversity gains to the data links, which 
are realized by using feedback CSI. Thus the CSI error causes the received power loss, which 
is quantified by the coefficients (1 — e^) and pi in Lemma [3l 

The distribution of the interference term |f^hmp of the SINR in Q is given in Lemma |4] that 
follows from [5, Lemma 2]. 

Lemma 4 ([5]): The random variable T„ „t := |frlhmp for n 7^ m is the product of the x^i^) 
random variable (Z + Q) and an independent L — 2) random variableO where Z and Q are 
identical to those in Lemma [3l 

'a l3(c,a/b) random variable T satisfies < T < 1 and has the CDF Pr(r < t) = b ^""^(l - z'')''''^dz with a > 0, 
6 > and z > [5]. 
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As observed from Lemma IH the number of scheduled queues (or spatial streams) K has no 
effect on the distribution of each interference power component 7e|f^hmp, but determines the 
number of such components (see Q). 

Using Lemma |4] and Lemma IH we obtain the loss on the departure rate for scheduled 
queues due to limited feedback. Define the departure rate conditioned on K = k as d{k) : = 
Ft (S\HR£ > 9 \ K = k) where SINR is given in Q. The loss on the departure rate can be 
characterized by the ratio d{k)/d{k). This ratio should be kept close to one by providing 
sufficiently accurate CSI feedback. Thus we define 5k := I — d{k)/d{k) and 

(5 := max(5fc = 1 — min d{k)/d{k) (13) 

which characterizes the maximum loss on the departure rate due to limited feedback. Given 5, 
the required number of CSI feedback bits per mobile is shown in the following lemma. 

Lemma 5: Given 5, it is sufficient for each mobile to provide the following number of CSI 
feedback back bits B 

B{5) = -{L-l)\og^5 + K (14) 
where k := (L - 1) loga [L{1 + LO) (l + |)). 

Proof: See Appendix |El □ 

Let C represent the stability region for zero-forcing SDMA with limited feedback. Using 
Lemma [51 the main result of this section is obtained and shown in the following theorem. 
Theorem 2: Given < 5 < 1, i?((5) in (fT4l) is sufficient for bounding the stability region C 

as (1 - 5)C C C C C. 

As observed from (fT4l) . B increases rather slowly (logarithmically) with the inverse of 5. For 
small values of 5, i? ~ (L — 1) logg |, where B increases linearly with (L — 1). In addition, 
B increases with the ratio 6/ P. Note that a large ratio corresponds to high coding rate or low 
transmission power and vice versa. 

B. Queueing Delay for Limited Feedback 

As mentioned, limited feedback reduces the departure rate for each queue and thereby increases 
queueing delay due to limited feedback. In this section, we derive B for bounding the degradation 
of the system delay performance. We consider both the Poisson and general packet arrival 
processes. 

1 ) Poisson Arrival: Without loss of generality, we consider the queue with the index 1 and 
the Poisson arrival rate A. Given a control policy tt specified by the probability vector p^, the 
average departure rate for perfect CSIT is defined as ji := [/i(p7r)]i where /i(p7r) is given in 
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Lemma [U Moreover, let fi denote the counterpart of /i for limited feedback. For Poisson arrivals, 
the average delay W follows the PoUaczek-Khinchin formula [22] 

W = - — ^—i- (15) 
2(1 -A/A) 

where X denotes the service time. The second moment E[X^] for ARQ transmission is given 
by E[X2] = (2 - fi)/fi^ [22]. Thus it follows from ([H that 

Let W denote the average delay for perfect CSIT. To limit the effect of CSI inaccuracy, it is 
desirable to bound the ratio W /W < M where M > 1. Given M, the required number of 
feedback bits B is shown in the following proposition. 

Proposition 4: Consider Poisson packet-arrival processes. Define r := 1 — ^. To satisfy the 
constraint W/W < M with M > 1, it is sufficient to provide B{6^) CSI feedback bits per 
mobile, where B{-) follows Lemma [5] and 5+ is given as 



2 



4(1-^7)- 



(17) 

□ 



;i + r)2 

Proof: See Appendix IB 

The dependence of M on B is illustrated in the following example. 

Example 2: In Fig. |4l the M versus B curves are plotted, where B is computed from M using 
Proposition m for L = 3, 6^ = 3 and P = 12 dB. Note that fractional values for B in Fig. |4] should 
be rounded to integers for practical applications. As observed from Fig. H as i? decreases, M 
increases very rapidly because fi approaches A. In the regime of large B, M converges gradually 
to 1 with increasing B. 

Finally, for — oo and thus M ^ 1, the relation between M and B is simplified in the 
following corollary of Proposition HI 

Corollary 1: For M — > 1, it is sufficient to scale B with M as 

5 = (L - 1) log2 + iL-l) log, ^^^^ + + O ((1 - 1/Mf) (18) 

where k is given in Lemma [51 

Example 2 (continue): In Fig. HI the curve of M versus the asymptotic approximation B = (L— 
1) log2 + K, from (fTSl) is plotted in dashed line. As observed from Fig. H the approximation 
is accurate for i? > 15 bit, corresponding to M < 1.5. 



Huang and Lau: Stability and Delay of Zero-Forcing SDMA with Limited Feedback 



15 




Fig. 4. The delay ratio M versus the number of feedback bits per mobile B (solid line) in Proposition |4] for Poisson arrivals, 
L = 3, 9 — 3 and P = 12 dB. The asymptotic approximation for Af ^ 1 is plotted in dashed line. 

2) General Arrival: In this section, we consider the same settings as in the preceding section 
but with general arrival processes. Let V represent the delay for a particular packet and thus 
= EfV"] . The effect of limited feedback on the CCDF of V is analyzed. Let X and Y denote the 
inter-arrival time and the service time, respectively. Using Kingman's bound for G/G/1 queues, 
the CCDF of V is bounded as [23] 

Ft{V >t)< exp{-f*t) (19) 
where f* is the solution of the following equation 

E[exp(r(f-X))] = L (20) 



As observed from (1191) . the probability that V exceeds the threshold t decreases exponentially 
with t, and f* determines the decrease rate. Let r* represent the counterpart of f* for perfect 
CSIT. Define the lose factor a := 1 — fi/ fi. Note that a < 6 with 5 in (fT3l) . The relation between 
f* and r* for a small S is established in the following proposition. Its proof uses perturbation 
theory [24]. 

Lemma 6: For cr — > 0, the CCDF exponential factors f* and r* for respectively limited 
feedback and perfect CSIT is related as 



r*-/(Oa + 0(a2) (21) 



where 



Proof: See Appendix |Gl □ 
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For high-resolution CSI feedback (B oo), we show that Kingman's bounds for perfect CSIT 
and limited feedback are related by linear scaling with a factor converging to 1 as i? ^ oo. This 
result is given in the following proposition. 

Proposition 5: For limited feedback with 

B = -{L-l)\og2r] + K + 0{r]), r] ^ (23) 

where k is given in Lemma |5J the delay CCDF is bounded as 

Fi{V >t) < e-''*\l + r]). (24) 
Proposition[5]is proved by combining (fT9l) and Lemma[6]and setting r] := f{r*)at. By integrating 
both sides of (l24l) . the average delay = E[V] can be bounded as 

W < r/ ^ 0. (25) 

Similarly, W <l/r*. Therefore, for B given in (l23l) . CSI inaccuracy increases the delay upper- 
bound by the factor (1 + r/). 

V. Numerical Results 

In this section, we present numerical results and further discuss the effects of limited feedback 
on the stability region and queueing delay of the SDMA system. 

Fig. [5] displays {1 — 5) with 5 in (fT3l) for different combinations of transmission power P 
and the number of transmit antennas L. The number of feedback bits per mobile is fixed at 
B = 12 bit and the SINR threshold 6* = 3. It is observed from Fig. [5] that (1 — 5) is insensitive 
to the changes on P for P > 5 dB. This is reflected in the expression of B in (fT4l) where B 
is independent of P for P ^ oo. Next, from Fig. [51 {1 — 5) is found to decrease rapidly as 
L increases. To keep {1 — 5) constant, B should increase linearly with (L — 1) as shown in 
(fT4l) . Thus L is an important factor for determining the CSI feedback requirement. It is worthy 
to mention that the same feedback requirement is found in [5] for the SDMA system with the 
ergodic capacity as the performance metric. 

Fig. [6] shows the departure rate per queue d{k) for different numbers of scheduled queues K 
and L = {2, 3, 4, 5}. The solid and dashed curves correspond to perfect CSIT and limited CSI 
feedback, respectively. We set {1 — 5) = 90% and B by evaluating (fT4l) . The vertical bars in 
Fig. [6] specify 10% rate loss with respect to the departure rate for perfect CSIT. It is observed 
from Fig. [6] that the rate loss due to limited feedback is within 5 = 10%. This confirms the 
sufficiency of B in (fT4l) . Moreover, Fig. [6] shows that for fixed L, the rate loss due to limited 
feedback is maximized at L scheduled queues. 
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Fig. 5. The loss factor (1 — 5) for the stability region for different combinations of transmission power and number of antennas. 
The number of feedback bits per mobile B = 12 bit and the SINR threshold 6 = 3. 



1.1 




Number of scheduled queues 



Fig. 6. Departure rate per queue for different numbers of scheduled queues. The solid and dashed curves correspond to perfect 
and limited CSI feedback, respectively. The vertical bars specify 10% rate loss with respect to the case of perfect CSIT. The 
number of antennas L = {2, 3, 4, 5} and the SINR threshold 9 = 3. 
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Fig. 7. Average queue length versus arrival rate per queue for both limited and perfect CSI feedback. The numbers of transmit 
antennas and scheduled queue are both equal to 4; the transmission power is P = 12 dB; for limited feedback, the number of 
feedback bit per mobile is B = {8, 10, 12, 20}; the SINR threshold is e = 3. 

Fig. |7] displays the curves of average queue length versus arrival rate per queue. The simulation 
parameters are L = A, K = A and P = 12 dB. The curves for limited feedback are plotted in 
solid line and that for perfect CSIT in dashed line. For limited feedback, B = {8, 10, 12, 20} 
bit and the SINR threshold 9 = 3. For given average queue length, the gain in arrival rate for 
additional feedback bits is more significant for smaller B. For example, for the average queue 
length of 50 packets, the gains are 0.6 and 0.4 for increasing B from 8 to 10 bit and 10 to 
12 bit, respectively. With B = 20 bit, the average queue length (or average delay) for limited 
feedback is close to that for that for perfect CSIT. 

VI. Conclusion 

In this paper, we consider the SDMA system with a L-antenna base station and L single- 
antenna mobiles. Given perfect CSIT and equal power allocation over scheduled queues, we have 
proved that the stability region of the SDMA system is a convex polytope. For any arrival-rate 
vector in this region, queues are stabilized by the control policy that maximizes the departure- 
rate-weighted sum of queue lengths. In addition, the stability region with power control has been 
derived. The stability region for case of limited feedback has been shown to be equal to the 
perfect-CSIT counterpart scaled by a factor that is determined by the number of CSI feedback 
bits per mobile B. We have also investigated the effects of limited feedback on the system delay 
performance. For Poisson arrival processes, limited feedback increases the average delay by a 
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linear factor M. For M — > 1, the required B has been shown to be = O (— log2(l — '^/M)). 
For general arrival processes, the tail probability of the delay for a particular packet can be 
bounded by Kingman's bound for perfect CSIT scaled by a factor (1 + r]). For — > 0, the 
required B has been proved to be = 0(— loggr/). The above results constitute a set of 
guidelines for designing practical SDMA systems with limited feedback and for delay-sensitive 
applications. 

This work also opens several issues for future research. First, the joint queue- and-beamforming 
control approach can be extended to SDMA systems with a large number of mobiles. Designing 
these systems should account for the effects of multiuser diversity gains on stability, delay and 
CSI overhead along the lines of [9], [10], [20], [37]. Second, the simplified CSI quantization 
model can be replaced by the actual codebook-based quantization [11]. Thereby the stability and 
delay related criteria can be included in codebook optimization. Finally, alternative beamforming 
algorithms such as minimum-mean- square-error (MMSE) or nonlinear precoding can be studied 
in the current framework. 

Appendix 

A. Proof for Lemma |2] 

From Lemma [U dT]) and dH]), the stability region C can be written as 

u u ^^^^ 

Let vi denote the ith element of TZ. Without loss of generality, let ri be the all-zero element 
of TZ. By the definition of TZ, (|26l) can be rewritten as 



j:,pi=lO<a<l 



E,P£=10<a<l 



L 



Define p[ = 1 — a Yle=2Pf^ P'e ~ ^Pf^ for 2 < £ < 2^. From (|27]) 

^ = U (28) 

Since > V£ and ^^p^ = 1, (l28l) shows that C is closed and contains all convex combinations 
of the elements of IZ. The desired result follows. 
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B. Proof for Theorem [7] 

Let C represent the polytope whose vertices belong to the set ext C in (flOl) . The proof for 
C = C comprises two steps: first, C is shown to be convex; second, the set 7^\ext C is shown to 
be a subset of C. 

Proving the convexity of C is equivalent to showing that its vertices are all exposed [36]. 
First, since all vectors in C except for are vectors of positive components, is exposed by the 
supporting hyperplane [36] 1-Lq = {v G : (v, 1) = 0} where 1 represents the vector of all-1 
components. Next, consider an arbitrary nonzero vertex s of C. We can write s = d(m)v for 
some m G X in (|9l) and v G Vm- Expand ext C as ext C := {ui, U2, ■ ■ ■ , u^-} where u„ G and 
T = I ext C\. An arbitrary point w G C can be writhen as w = J2^=i '^nU„ where a„ > V n and 
J2n=i = 1- Thus (w, v) = Yln=i '^ri(u„, v). The term (u.„, v) belongs to one of the following 
four cases: 

1) If u„ = (i(m)v, then (u„, v) = md{m) since v G Vm', 

2) If u„ G TZm and u„ 7^ d{m)\, then (u„, v) < ■md{m); 

3) If u„ G TZa with a > m, then (u„, v) < md{a) < md{m); 

4) If u„ G TZa with a < m, then (u„, v) < ad{a) < ■md{m) where the last inequality follows 
from dH). 

By combining these cases, we conclude that (u„, v) < md{m) and thus {w, v) < md{m) if 
w 7^ s. In other words, C lies in a half-space determined by the hyperplane Tig = {a G : 
(a, v) = m(i(m)}, and Tig^C = s. Thus s is exposed by Tig- Since s is arbitrary, it follows that 
all vertices of C are exposed, proving the convexity of C. 
Define the complementary set of X in Q as 

J= = J < A; < L : kdik) < max mdim) \ . (29) 

1^ 0<m<k J 

Note that 7^\ext C = Ufcej^ ^fc- Given that both C and C are convex polytopes, to prove C = C, 
it is sufficient to show that the set [jj^^^jc Tlk ^ C. Specifically, we show that an arbitrary point 
d G UfcGX^ belongs to C as follows. There exists a E such that d G TZa and d = (i(a)v 
where v G V^. Define m := arg maxo<m<a . From (l29l) . ad{a) < md{m). Furthermore, 
from (0), fh el and thus TZfh G ext C. Next, define the set W := {v G V™ : (v, v) = m}. Note 
that W is nonempty as v G Va and v G Vm and a > m. Moreover, define z := Zlvew d{rh)v 
that is the convex combination of the points in TZm, which are the vertices of C. Since C is a 
convex polytope, z G C. It follows that the line segment [0, z] G C. Note that z = ^d{m)v. 
Since d{a) < ^d(m), d G [0, z] and thus d G C. Thus IJfcex^ ^ ^- "^^^^ completes the proof. 
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C. Proof for Proposition [7] 

Let dk{P) denote the same departure-rate function as d^iy) except with the different argument 
P. with the Consider k> 1 and m < k. Using Lemma [T] and L'Hospital's rule 

lim \ = lim 
= lim 



r(L- 


- m + 1) /m ^e" 


-"da 


r(L 


-k + 1) jZ h^-^e- 


Mb 


r(L- 


- m + 1) \'m a^~^e~ 


-"da 


r(L 


- A; + 1) jZ h^-^e- 


-Mb 


r(L 


-m+1) 


kB 
P 



P^o T{L-k + 1) (m^)^-™ 
= 0. (30) 

It follows that the index set X defined in ^ converges with P ^ oo as 

limX(P) = 7^l U7^o = 7^l U {0}. (31) 

The first claim in the proposition statement follows from Theorem [T] and (|3TI) . 
Using Lemma [Hand Alzer's inequalities [38], d^ can be bounded as 

<dk<l- il~e-—) (32) 



where f3 = [r(L - A; + 1)] ^-k . For m< k and using ([32|), for P > 1 



4(P) - jdmiP) 



L-k 





2{L-ky 

-+011^1 1,(33) 

From ([33l) and ©, limp^oo^(-P) = Combining this result and Theorem [T] completes the 
proof for the second claim in the proposition statement. 

D. Proof for Lemma \3\ 

The result for K = L follows from the fact that f„ and h„ are independent [5]. The proof for 
the case oi K < L h given as follows. Without loss of generality, consider the random variable 
Let Ti denote the {K — l)-dimension vector sub-space containing (h2,h3,--- i^.r)- 
A basis of H comprises the vectors (ei,e2, ■ ■ ■ ,ex-i). The null space of H, represented by 
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Ti^, is spanned by the basis vectors {eK^GK+i^^h)- Given these definitions, the unitary vector 
Si := hi /I hi I can be decomposed as Si = Xl^^i -^i/e^ where Yld=i l-^i/P = 1- For zero-forcing 
beamforming discussed in Section III-BI 

2 





2 


L 


2 


L 


f^Si 


= max 


1=1 


= max 





It follows that 



fi 



(34) 



Define Si := hi/||hi|| and decompose it as Si = X]t=i -^i/^^- The projection of Si into the 
subspace is hence Si := J2e=K+i -^i/^e- Moreover, decompose hi as hi = J2e=i 
Similarly, the projection of hi into is given as hi = J2m=K+i ^hm^rn- Using above definitions 
and dMl) 



hlfi 



■51, 



H-^ ,2 UU l|2rt-^ |2 

|SiSi| > ||hi|| |s|si| . 



(35) 



The inequality (a) follows from J2£=k I'^i/P < 1- Based on quantized CSI model in Section HTCl 
and by geometry [39] 



i£ta |2 

SiSi 



|hi| 



|hi| 



Pi- 



|hi| 



|hi| 



|hi| 



where — ei < pi < ei. From (1351) and (136 



hlfi 



>(l-ei)||hif + p(||hi 



|hil 



(36) 



(37) 



Given that {hi m} are i.i.d. exponential random variables with unit mean, the desired result 
follows from (l37l). 



E. Proof for Lemma |5] 

Let A denote the set of indices of the K scheduled mobiles and nonempty queues. To simplify 
notation, define Ts = Ylim&AT^ m- Since {T^™} are i.i.d. P{l,L — 1) random variables from 
Lemma S Ts < - 1 and E[Ts] = 1. Using Lemma S Lemma H and (g]), the SINK can be 



rewritten as 

SINR. 



1 + -feeiZe + Q,)Ts 



- l + 7e,Ts(Z, + Q,)' 
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From (l38l) . the departure rate d{K) is lower bounded as 

V 1 + 'jTY;ee{Ze + Qe) 

9h + e,Qfil + T^> 



E[Pr(Z > 



= EK:o"n"exp(-n) 

> E[5^^'f(^/7)"exp(-n)" 



E 



L-K 



> E 



> 



[exp (-n + eh)] J2 \ (o/ire 

eeil+T^e){Qe + e/^) \ 



exp 



l-e,(l + Ts^^) ; 
E[6,](l + Tsg)(E[g,]+g/7) 

1-2^^(1 + (is: - 1)^) 



d{K) 
d{K) 



1_2-— (l + (ir- 1)0) 



f{B) 

The equality (a) holds since Zi is a x(L — K + 1) random variable according to Lemma HI The 

B 

inequality (b) is obtained using eg < 2 ^-i and < K — 1. The equality (c) follows from 

E[ei] = 2^T^ (1 - i)' E[Ts] = 1 and E[Qe] =K -l.By setting f{B) < 6, we obtain that 

B ^ 6 

- {l-i)il + e){K~l + ^)+6il + {K- 1)9) ■ 

Let K = L and we have more stringent feedback requirements 

B ^ 6 

- (L-1)(1 + ^^)(1-1 + |)+5(1 + (L-1)0)- 

The desired result follows from the above inequality. 



F. Proof for Proposition |?] 

The average queueing delay W is given by (fT6l) but with /t replaced by fi. From (fT6l) and ([13 

1 



W W 

14? ~ l-5"/t-A 



u — A 2 — u 

X X 

2-/i 



< 



1 



X 



/i — A 

1 — 5 /i — A 

r 

:i-s){T-6)- 



(39) 



Huang and Lau: Stability and Delay of Zero-Forcing SDMA with Limited Feedback 



24 



where (a) holds since fi/ fi > 1 — 5 from (fT3l) and the definitions of jj, and (i. Therefore a sufficient 
condition for W/W < M h 



^ < M 



(l-5)(r-5) 



or equivalently 



It follows that 



5^-(l + r)5+(l-— )r>0. (40) 



^ ^ l + ^-Va+^)^-4(l-iT)r 



1 + r 



(1 + r)^ 



(41) 



n=l 

r 

(42) 



One can check that the right-hand- side of (|4T)) is always real for M > 1. The desired result 
follows from ((4T)) and Lemma [5l 

G. Proof for Lemma |6| 

Given ARQ transmission, the moment generating function of Y is obtained as 

oo 

E[e''^] = ^/i(l-/i)"-ie™ 

Since F and X are independent, by substituting (|42l) into (|20|) 

/iE[e-^*^] - e""* + 1 - /i = 0. (43) 

We can write /t = (1 — cr)/i with < cr < 1. Note that a < 5 with 5 given in (fT3l) . Moreover, 
define r* as the solution of the following perfect-CSIT counterpart of (143]) 

/xE[e-"*^] - e""* + 1-/2 = 0. (44) 

For 5^0 and hence a ^ 0, the asymptotic relationship between f* and r* is obtained using 
perturbation theory [24] as follows. Expand f* as f* := a + [3a + 0{a'^). By substituting this 
expression into (|43l) 

/iE[e-"^] - e-° + 1 - - {/iE[e-"^X]/? + /iEfe"'"^] - /Je"" - /i} (x + 0{a'^) = (45) 



Hi n2 
Following the approach of perturbation theory, the terms Hi and 112 are set equal to zero, thus 

/iE[e-"^] - + 1 - = 

/xE[e""^X]/? + /iE[e-"^] - /^e"" - = 0. 
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By solving the above equations, we obtain that a = and P := —f{r*) with /(r*) in (|22l) . 
This completes the proof. 
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